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Abstract 

Solutions of a variational inequality are found by giving conditions 
for the monotone convergence with respect to a cone of the Picard 
iteration corresponding to its natural map. One of these conditions 
is the isotonicity of the projection onto the closed convex set in the 
definition of the variational inequality. If the closed convex set is a 
cylinder and the cone is an extented Lorentz cone, then this condition 
can be dropped because it is automatically satisfied. The obtained 
result is further particularized for unbounded box constrained varia¬ 
tional inequalities. For this case a numerical example is presented. 
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1 Introduction 


In this paper we will study the solvability of variational inequalities on closed 
convex sets by using the isotonicity of the metric projection mapping onto 
these sets with respect to the partial order defined by the cone. Apparently 
this approach has not been considered before. 

Variational inequalities are models of various important problems in physics, 
engineering, economics and other sciences. The classical Nash equilibrium 
concept can also be reformulated by using variational inequalities. 

It is known that a vector is a solution of a variational inequality if and only 
if it is a zero of the corresponding natural mapping pQ. Hence, a variational 
inequality on a closed convex set is equivalent to finding the fixed points of the 
difference between the identity mapping and its natural mapping. The latter 
mapping is just the composition between the projection onto the latter closed 
convex set and the difference between the identity mapping and the mapping 
dehning the variational inequality. I we could guarantee the isotonicity of the 
mappings in the latter composition with respect to the partial order defined 
by the cone, then by an iterative process we could construct an increasing 
sequence with respect to the partial order defined by the cone by using Pi¬ 
card’s iteration. If we could also guarantee that this sequence is bounded 
from above with respect to the partial order defined by the cone, then this 
sequence would be convergent to the solution of the variational inequality (in 
fact it would be convergent to a solution of the corresponding equivalent fixed 
point problem). It turns out that there is a class of variational inequalities 
and a class of cones, that extend the Lorentz cone, for which this idea works 
very well. The only restriction is that the variational inequality has to be de¬ 
fined on a cylinder. Such problems appears in the practice. For example the 
unbounded box constrained variational inequalities are of this form. Based 
on the above idea a theorem for finding solutions of variational inequalities 
on a cylinder will be presented and an example will be given. Similar ideas 
to the above ones were presented in [Ml for complementarity and implicit 
complementarity problems, but with a strong restriction on the cone defining 
the problem. The idea of monotone convergence for complementarity prob¬ 
lems defined by general cones is considered first in the paper [8] . The present 
paper extends the results of [8] for variational inequalities. 

Several other papers dealt with conditions of convergence for iterations 
similar to the above one, as for example P-HU]. However, neither of these 
works used the ordering defined by a cone for showing the convergence of the 
corresponding iterative scheme. Instead, they used as a tool the Banach fixed 
point theorem and assumed Kachurovskii-Minty-Browder type monotonicity 
(see [201123] I and global Lipschitz properties of F. 
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The structure of the paper is as follows: In the Preliminaries we intro¬ 
duce the terminology and notations use throughout this paper. In Section 
3 we recall the definition and basic properties of the extended Lorentz cone, 
and determine all sets onto which the projection is isotone with respect to 
this cone. In Section 4 we will find solutions of a variational inequality by 
analyzing the monotone convergence with respect to a cone of the Picard 
iteration corresponding to its natural map. In Section 5 we will particular¬ 
ize these results to variational inequalities dehned on cylinders, by using the 
extended Lorentz cone for the corresponding monotone convergence above. 
In this case we can drop the condition of Proposition [T] that the projection 
onto the closed convex set in the definition of the variational inequality is 
isotone with respect to the extended Lorentz cone, because this condition is 
automatically satisfied obtaining the more explicit result of Theorem [21 The 
latter result extend our results in [8] for mixed complementarity problems. In 
Section 6 we further particularize the results for unbounded box constrained 
variational inequalities. In Section 7 we will give a numerical example for 
this case. 


2 Preliminaries 

Denote by N the set of nonnegative integers. Let k,m,p,q E N\ {0} and 
be the u-dimensional real Euclidean vector space. 

Dehne the direct product space MP x as the pair of vectors {x, u), where 
a; G R^ and u G R^. 

Identify the vectors of R”^ by column vectors and consider the canonical 
inner product in R™' defined by R”^ x R™' 3 {x, y) ^ x~^y G R with induced 
norm R™ 3 a: ||a;|| = \Jx~^x G R. 

The vector space R^ x R'^ can be identified with via the mapping 
RP X R'^ 9 {x,u) I—)■ {x'^,u'^)~^ G RP+'^. This identihcation leads to the inner 
product in R^ x R"? given by 

((a;, u), {y, v)) x~^y -|- u^v 

The closed set K C R'" will be called a cone if 77 fl (—77) = {0} and 
Xx + yy E K, whenever X, y > 0 and x,y E K. Let 77 C R^+^ be a cone. 
Denote <k the relation defined by x <k y y — x E K and call 

it the partial order defined by 77. The relation <x is reflexive, transitive, 
antisymmetric and compatible with the linear structure of R'" in the sense 
that X <K y implies that tx + z <k ty -|- z, for any z G R™' and any f > 0. 
Moreover, <k is continuous at 0 in the sense that if a;"' —)■ a; when u —)■ oo and 
0 <K for any u G N, then 0 <k x. Conversely any reflexive, transitive and 
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antisymmetric relation ^ in R"* which is compatible with the linear strnctnre 
of R™ and it is continnons at 0 is dehned by a cone. More specihcally, :<=<k, 
where K = {x ^ R™ : 0 R a:} is a cone. 

For any closed convex set C denote by Pc '■ R™ —t R™ the metric projec¬ 
tion mapping onto C, that is, the mapping dehned by Pc{x) G C and 

||x - Pc{x)\\ = min{||a: - y\\ ■. y & C}, 

for any x G R™'. Since C is closed and convex, the projection mapping is 
well dehned and by its dehnition it follows that 

Py+c{.x) = y + Pcix - y) (1) 

for any x, ?/ G R™. It can be also shown that Pc is nonexpansive (see m, 
that is, 

\\Pc{x) - Pcm < \\x-yl (2) 

for any x, y G R”^. 

Let K C R™ be a cone. The mapping F : R”* —)■ R™ is called K-isotone 
if X <K y implies F{x) <k F{y)- 

The nonempty closed convex set C C R"* is called iF-isotone projection 
set if Pc is iF-isotone. 

The set C R™ is called iF-bonnded from below [K-bounded from above) 
if there exists a vector y G R™ snch that y <k x (x <k y), for all x G 12. In 
this case y is called a lower K-bound {upper K-bound) of 12. If y G 12, then 
y is called the K-least element {K-greatest element) of 12. 

Let X C N be an nnbonnded set of nonnegative integers. The seqnence 
{x^}n£i is called iF-increasing (LF-decreasing) if x^^ <k x"'^ (x"'^ <k x"'^), 
whenever rii < n 2 - 

The seqnence is called K-bounded from below {K-bounded from 

above) if the set {x^ : n G X} is iF-bonnded from below (iF-bonnded from 
above). 

A closed convex cone K is called regular if any iF-increasing seqnence 
which is iF-bonnded from above is convergent. It is easy to show that this 
is eqnivalent to the convergence of any iF-decreasing seqnence which is iF- 
bonnded from below. It is known (see [25]) that any cone in R™ is regnlar. 

The dual of a cone K C R”^ is the cone K* dehned by 

iF* = {x G R”^ : x~^y = 0, V|/ G iF}. 

A cone K C R™" is called polyhedral if it is generated by a hnite nnmber 
of vectors v^,... ,v^, that is, 

K = cone{n\ ..., := {AiX^ H-h : Ai,..., A^ > 0}. 
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The vectors v^,... ,v^ are called the generators of K. 

The affine hyperplane with normal u G \ {0} and through a G is 
the set defined by 

T-L{u, a) = {a: G : {x — a,u) = 0}. (3) 

An affine hyperplane 'H(u,a) determines two closed halfspaces 'H-{a,u) and 


l-i+fa^ a) 

of R™, defined by 




l-L-{u^ a) 

= {x G R’" : 

(x — a,u) < 0}, 

and 

TL+iu, a) 

= {x G R™ : 

(x — a,u) > 0}. 


3 Extended Lorentz cones 

Let p, q be positive integers. For a, 6 G denote a > 6 if and only if b <rp a, 
that is, components of a are at least as large as the corresponding components 
of b. Denote by e the vector in with all components equal to one and by 
e* the canonical unit vectors of In [8] we defined the following notion of 
an extended Lorentz cone: 

L = {{x,u) G X : X > ||M|le} (4) 

and showed that the dual of L is 

L* = {{x,u) G R^ X R^ : x'^e > ||m||,x > 0}. 

In the same paper we also showed the followings: 

• The extended Lorentz cone L defined by (jl]) is a (regular) cone. 

• The cone L (or L- 123) is a polyhedral cone if and only if g = 1. 

• If g = 1, then the minimal number of generators of L is 

{p + 2)(1 — (5pi) + 2(5pi, 

where 6 denotes the Kronecker symbol. 

• Ifg = l,p = l, then a minimal set of generators of L is 

{(1,1), (1,-1)}, 

and if g = 1, p > 1, then a minimal set of generators of L is 
{(e, l),(e,-l),(e*,0) : i = l,...,p}. 
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• If g = 1, then L* is a p + 1 dimensional polyhedral cone with the 
minimal number of generators 2p and a minimal set of generators of L* 

is 

{(e*,l),(e*,-l) : i = 

• If g = 1 and p > 1, then note that the number of generators of L and 
L* coincide if and only if they are 2 or 3-dimensional cones. 

• The cone L is a subdual cone and L is self-dual if and only if p = 1, 
that is, L is the g -|- 1-dimensional Lorentz cone. 

• L is a self-dual polyhedral cone if and only if p = g = 1. 

In Theorem 2 of [S] we determined the L-isotone projection sets. For 
convenience we repeat this theorem here: 

Theorem 1. 

1. Let K = ML X C, where C is an arbitrary nonempty closed convex set 
in M*? and L he the extended Lorentz cone defined by (SI)- Then, K is 
an L-isotone projection set. 

2. Let p = 1, g > 1 and K G M^ x M'^ be a nonempty closed convex set. 
Then, K is an L-isotone projection set if and only if K = M^ x C, for 
some C gM'^ nonempty closed convex set. 

3. Let p, g > 1, and 

K = C X Mfi 

where Y = is a unit vector. Then, K is an L-isotone projection 

set if and only if for each i one of the following conditions hold: 

(a) The vector = 0. 

(b) The vector = 0, and there exists i j such that af = \f2j2, 
aj = —\/2/2 and al = 0, for any k ^ {*, j}- 

4 Variational inequalities 

Let K G M™' be a closed convex set and F : —)■ be a mapping. Then, 

the variational inequality VI{K,F) dehned by F and K is the problem of 
hnding an x* G iF such that for any y E K, 

{y — xYF{x) > 0. 
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It is known that x* is a solntion of VI{F, K) if and only if it is a fixed point 
of the mapping I — = Px o [I — F), where I is the identity mapping 

of MF and is the natnral mapping associated to VI{F,K) defined by 
Fx^^ = I — Px o (/ — F) [1] . Consider the Picard iteration 

= Px{x'^ - F{x^), (5) 

If F is continnons and is convergent to x*, then by a simple limiting 

process in ([5j), it follows that x* is a fixed point of the mapping Px o (/ — F) 
and hence a solntion of ld/(F, K). Therefore, it is natural to seek convergence 
conditions for x". Let us first state the following simple lemma: 

Lemma 1. Let K C he a closed convex set, F : MF —)■ R™' be a continu¬ 
ous mapping and L be a cone. Consider the sequence {a:"'}neN defined by (jS]). 
Suppose that the mappings Pk and I — F are L-isotone, x^ <l x^ , and there 
exists a y E R™' such that x"' <l y, for all n E sufficiently large. Then, 
is convergent and its limit x* is a solution of VI{F, K). 

Proof. Since the mappings Pk and I — F are L-isotone, the mapping x i—)■ 
Pk o (J — F) is also L-isotone. Then, by using ([5]) and a simple inductive 
argument, it follows that is L-increasing. Since any cone in R™' 

is regular, {x"'}„gi^ is convergent and hence its limit x* is hxed point of 
Px o (J — F) and therefore a solution of VI{F, K). □ 

Remark 1. The condition x° <l x^ in Lemma [1] is satisfied when x^ E 
KnF-\-L). Indeed, z/x° G KnP-fi-L), then -F{x^) E L and x^ E K. 
Thus x° <L x° — F(x°), and hence by the isotonicity of Px we obtain x° = 
Pk{,x^) <l Pk{,x^ - F{x^)) = x\ 

Proposition 1. Let K C R”^ he a closed convex set, F : R™ —)■ R™ he a 
continuous mapping and L be a cone. Consider the sequence {x”}„gN defined 
by dS]). Suppose that the mappings Px and I — F are L-isotone and x° <l x^ . 
Denote by I the identity mapping. Let 

= {x e F n (x° + L) : F(x) G L} 

T = {x G F n (x° -I- L) : Px{x — F(x)) <l x}. 

Consider the following assertions: 

(i) 

(ii) T ^ 0, 
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(in) The sequence {x^'jneN is convergent and its limit x* is a solution of 
VI{F,K). Moreover, x* is the L-least element ofT. 

Then, C F and ^ ^ 

Proof. Let us first prove that C F. Indeed, let y E fl. Since Pk is L- 
isotone, y — F{y) <l U implies PKiy - F{y)) <L PRiy) = U) which shows 
that y eT. Hence, H C F. Thus, ([I]) (jll]) is trivial now. 

© ^ dmD; 

Suppose that F 7 ^ 0 . Since the mappings Pk and I — F are L-isotone, the 
mapping Pk o [I — F) is also L-isotone. Similarly to the proof of Lemma [H 
it can be shown that is L-increasing. Let ?/ G F be arbitrary but 

hxed. We have y — x^ E L, that is x^ <l y. Now, suppose that x” <l y. 
Since the mapping Pk o [I — F) is L-isotone, x” <l y implies that x^^^ = 
Pk^x"^ — F{x"‘)) <L Pxiy — F{y)) <l y. Thus, we have by induction that 
x'^ <L y for all n G N. Then, Lemma implies that is convergent and 

its limit X* E Kr[ (x° + L) is a solution of VI{F, K). Since x* is a solution of 
VI{F,K), we have that Pk{x* — F{x*)) = x* and hence x* E F. Moreover, 
the relation x^ <l y in limit gives x* < y. Therefore, x* is the smallest 
element of F with respect to the partial order dehned by L. □ 


5 Variational Inequality on cylinders 

Let p, q be positive integers and m = p + q. By a cylinder we mean a set 
K = M.P X C <Z W xM'^ = M™. In this section we will particularize the results 
of the previous section for variational inequalities on cylinders. 

Lemma 2. Let K = x C, where C is an arbitrary nonempty closed convex 
set in M^. Let G : W x W ^ H : x —)■ M”? and F = {G,H) : 
Rr’ X —>■ X R'^. Then, the variational inequality VI{F, K) is equivalent 
to the problem of finding a vector (x,u) G R^ x C such that 

I G(x,u) = 0, 

[ (v — u)~’~H(x,u) > 0 


for any v E C . 

Proof. The variational inequality VI{F, K) is equivalent to hnding an [x, u) E 
Rr’ X C such that 


{y — xyG{x, u) + {v — H{x, u) > 0 


( 7 ) 


for any {y^v) e MP x C. Let {x,u) e x C be a solution of ([7]). If we 
choose V = u E C in ([7j), then we get {y — x)~^G{x,u) > 0 for any y G 
Hence, G{x,u) = 0 and (n — u)~^H{x,u) > 0. Conversely, if {x,u) G x C 
is a solution of ([6]), then it is easy to see that it is a solution of ([7j). □ 

By using the notations of Lemma (|2]) the Picard iteration ([S]) can be 
rewritten as: 

r ^n+l ^ x^-G{x^,u^), 

I ^n+l ^ Pc[u'^ — H{x‘^,u"')). 

Consider the partial order defined by the extended Lorentz cone (jlj). Then, 
we obtain the following proposition. 

Theorem 2. Let K = x C, where G is a closed convex set. Let G : 
M^xM'? —)■ H : RP X —)■ be continuous mappings, F = {G, H) : M^x 
RP X R^. Let G R^ x C and consider the sequence ( 2 :", 

defined by Letx,y andu,v . Suppose that x^—x^ > \\u^—u^\\e 
(in particular, by Remark{l\ this holds if —G{x^,u^) > \\H{x^,u^)\\e) and 
that y — X > ||n — M||e implies 

y — X — G{y,v) + G{x,u) > ||n — u — H{y,v) + H{x,u)\\e. 


Let 

n = {{x,u) eWxG:x — x^> ||m — M°|le , G{x,u) —x^> \\H{x,u) —u^\\e} 
and 


T = {(x, u) G R^ X C : X — > ||m — u^\\e, 

G{x, m) — x° > ||m — — Pc{u — H{x, M))||e} 

Consider the following assertions 

(i) 

(ii) T ^ 0, 

(Hi) The sequence {(x ”,is convergent and its limit {x*,u*) is a so¬ 
lution ofVI{F,K). Moreover, {x*,u*) is the smallest element ofT with 
respect to the partial order defined by L. 

Then, hi C T and ^ ^ 
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Proof. Let L be the extended Lorentz cone defined by (jl]). First observe that 
K n (x° + L) 7 ^ 0. By nsing the definition of the extended Lorentz cone, it 
is easy to verify that 

n = Kn ((x°, + L) n = {zeKn ((x°, + l) ■. f{z) g l} 

and 

T = {zEKn ((x°, n°) +L)-.Pk{z- F{z)) <l z}. 

Let x,y eW and u,v E C. Since y — x > ||n — M||e implies y — x — G{y, v) + 
G{x,u) > \\v — u — H{y,v) + H{x,u)\\e, it follows that I — F is L-isotone. 
Hence, by Proposition [T] (with m = p + q) and Lemma [21 it follows that 
C r and ([i]) ([n]) dm]). □ 

6 Unbounded box constrained variational in¬ 
equalities 

Let p, q be positive integers, m = p + q and K = bi] be a box, where 

ai, 6^ G M U {—oo, cxd} and < bi, for all £ G {1,..., m}. The f-th entry of 
the projection fnnction is (see Example 1.5.10 in my- 


{ Qji if Xj ^ flj, 

Xi if ai<Xi< bi, (9) 

bi if bi < Xi. 

Let B = y.^^i[ai,bi\ C and G = 6 p+j]. So we have 

PK{y^v) = {PB{y),Pc{v)) ( 10 ) 

and the Picard iteration (|S|) becomes 

x”+^ = mid(ai, h, (x*" - F(x”))i). (11) 

Let L be the extended Lorentz cone defined by (jlD- The next proposition 
shows that the L-isotonicity of a box is eqnivalent to the box being a cylinder. 

Proposition 2. Let L he the extended Lorentz cone defined by (jlD- Then, 
the projection mapping P^ is L-isotone if and only if K = x G where 
G = [®p+j) ^p+j] ■ 
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Proof. The sufficiency follows easily from item 1 of Theorem 2 in [S] (repeated 
in the Preliminaries as Theorem [T]) For the sake of completeness we provide a 
proof here. Suppose that B = W. If (x, u) <l {y,v), that is, y—x > ||n—M||e, 
then by the nonexpansivity of Pc we get 

Psiy) - Pb{x) = y - X >\\v - u\\e > ||Pc(^^) - Pc{u)\\e, 
which is equivalent to 

Pk{x,u) = {Pb{x),Pc{u)) <l {PB{y),Pc{y)) = PK{y,v). 

Hence, Pk is L-isotone. Although, the necessity could also be derived from 
item 3 of the same theorem, it is more clear to prove this directly as follows. 
Suppose that Pk is L-isotone. We need to prove that Oi = —cx), hi = oo, 
for any i = 1,.. .p. Assume to the contrary, that there exist at least one 
k G {1,... ,p} such that either or bk is a hnite real number. Assume that 
bk is a hnite real number. The case is a hnite real number can be treated 
similarly. Let u and v be two diherent vectors in C. Then, Pc{u) = u and 
Pciv) = V. It is easy to choose x,p G such that yi — Xi > ||n — m|| for 
all i G {1,... ,p} and bk < Xk < yk- For example, we may choose Xj = 6ikbk 
and yi = Sikbk + ||n — m||, for all i = {1,... ,p}, where 6ik is the Kronecker 
symbol. Then, (x,m) <l {y,v) and by Q we have {PK{y,v))k = (Pi^(x,n))fc, 
or equivalently {PB{y))k = (^^(x))^. Hence, by flTnll and the L-isotonicity 
of Pk we get 

0 = (Pb(p) - PBix))k > \\Pciv) - Pciu)\\ = 11^^ - u\\ > 0, 

which is a contradiction. 

Hence, the results of Theorem [2] can be particularized to the set K given 
by Proposition [21 with the Picard iteration taking the form fllll) and the 
function ’mid’ given by ([9]). In the next section we will present an example 
for this particularized result. □ 

7 Numerical example 

Let K = X C where C = [—10,10] x [—10,10]. Let L be the extended 
Lorentz cone dehned by dl]). Let/i(x,-u) = 1/12(xi + ||m||+ 12) and/ 2 (x, u) = 
1/12 (x2 + ||m|| — 7.2). Then it is easy to show that these two functions are 
L-monotone. Let wi = (1,1,1/6,1/3) and W 2 = (1,1,1/3,1/6) so wi and W 2 
is in L. For any two vectors {x,u) and {y,v) in K, suppose (x,m) <l {y,v), 
we have Pi — Xi > ||n — m|| > IlnH — ||n|| by triangle inequality. Hence, 

fiiy,v) - /i(x,m) = -xi - (||m|| - ||n||)) > 0. 
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Similarly, we can prove that if {x,u) <l (|/,n), then f2{y,v) — f2{x,u) > 0. 
Provided that K is convex, and Wi,W2 G L, if {x,u) <l {y,v) holds, then 

ifiiy^v) - fi{x,u))wi + if2{y,v) - f2ix,u))w2 G L. 

Thns, fi{x,u)wi + f2{x,u)w2 <l fi{y,v)wi + f2{y,v)w2. Therefore, the 
mapping fiWi + /2W2 is T-isotone. Hence, choose the fnnction 

{X — G,U — H) = fiWi + /2W2 = ^/i + / 2 , fi + / 2 , -/i + -/ 2 , -fi + g, 

( 12 ) 

where G, H, fi and /2 are considered in the point {x,u). It is necessary to 
check that all the conditions in Theorem [2] are satished. First, let {x^,u^) = 
(43/30,13/30, 2, 5) which is in iP = x G. Consider (x, u) = (15,15, 6 , 8 ). 
Thns, 


X — X 


0 


/407 437\ 

V"3(r’^ ) 


> (5,5) 


11(4,3)11(1,1) = ||n-n°||e. 


By the dehnition of the mapping {G, H), we know that G{x, u) — xo = x — 
ifi + /2)e - x° and H{x,u) - = u - (I/6/1 + I/3/2, l/ 3 /i + I/6/2) - 

Then, 

1 /43 13\ 

Gix, u)-xo = (15,15) - -(15 + 15 + 2 X 10 + 4.8)e - (- j = ( 9 ,10) 


and 

/121 91 \ 

7/(x,n)-n° = (6,8)-(^—,-j-(2,5) 
Hence, we can easily see that: 


359 91\ 
40 ) 


G{x, u) — Xq > II (3, 3)11 e > ||i7(a;, n) — M°||e, 


which shows that (x, n) G and 7 ^ 0 . Therefore, the conclnsions of 
Theorem [2] will apply for onr example. 

Now, we begin to solve the VI. Snppose that (x, n) is its solntion. Since 
G(x, u) = 0 , and 

X - G{x, u) = (/i + / 2 , /l + / 2 ), 

where fi = fi{x, n), i = 1, 2, we have Xi = X 2 = /i + / 2 - Moreover, since 

xi = ^(a:i + X 2 ) + ^||m|| + 0.4, 

12 0 


12 



we get 


( 13 ) 


1 „ „ 12 

X. = .. = -Hull + 


Obviously, if (x,u) G ker(F), that is, G{x,u) = 0 and H{x,u) = 0, then 
{x, u) is a solution of the variational inequality. The equality H{x, u) = 0 
implies 

Ul = |/i + I/ 2 , 


U2 = Ifl + 1 / 2 - 

By using equations flT^ . equations ffTTl) become 


Ml = 


20 


\u\ 


75’ 


“2 = ^I|m|| + |I- 


Hence, U 2 = ui 


20 II “'ll ' 75' 

By substituting back to ffT5|) . we obtain 
28 


Ml Ml 


995 


= 0 . 


However, the equation ffT5|) can be transformed to 

19mi + M2 = ||m|| 

By squaring both sides, we get 


(14) 


(15) 


361mi + M2 + 38miM2 = u\ + u\. 


Hence, 


2mi(180mi + I9M2) = 0 


(16) 


Thus, if Ml = ^ , then since M 2 = Mi + ^ equation (ITB)) will not hold. 
Therefore, the only solution in this case is 


(a;,M) 


—, 0 ,— 

15’15’ ’15 


Next, suppose that the solution (x, m) ^ ker(F), that is, H = {x,u) 7^ 0. 
Then, by (E]) and flTSD . we get 

(mi - mi)(300mi - 15||m|| + 4) + (m2 - M2)(300m2 - 15||m|| - 76) > 0 (17) 
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From the general theory of variational ineqnalities, it is known that u 
shonld be on the bonndary of C. Let Ui = 10. By choosing vi < Ui and 
V2 = U2 and nsing ||m|| < 20, it can be seen that ineqnality flTTD will not 
hold. Similarly, if U 2 = 10, by choosing V 2 < U 2 and vi = mi, ineqnality f[T7|l 
will not be satished. If u\ = —10, then by choosing vi > ui and V 2 = M 2 , 
ineqnality 0171) will not hold. Similarly M 2 = ~10 will not lead to a solntion. 
Hence, this case will not provide a new solntion. Therefore, the only solntion 
of the variational ineqnality is 

obtained above. 

The Picard iteration can be completed by nsing any spreadsheet software. 
Note that since the variational ineqnality is box constrained, the iteration 
in (ITT|) will be calcnlated by nsing the median fnnction as it is shown in the 
following tables. More precisely, the initial point is given in the hrst row 
for n = 0 and the later iterations are given in the following rows, where the 
colnmns below n" and obtained by the median of the npper bound, 

lowerbound and u'^ — H{x"',u‘^). In the hrst table, we showed the iteration 
using the previous example where = (43/30,13/30,2,5). In the 

other tables, we showed the iteration from different initial points in different 
directions outside of C. It can be observed that the iterations in all tables 
converge to the same unique solution of the variational inequality. 


n 

Xi 

rf^Tt 

^2 


U2 

0 

colo 

13 

30 

2 

5 

1 

1329 

1329 

244 

15 

865 

865 

973 

29 

2 

3 

494 

10 

9 

4 

657 

183 

28 

3 

40 

575 

5 

7 

69 

992 

432 

26 

4 

19 

284 

0 

4 


35 

523 

15 

5 

53 

53 

0 

4 


99 

99 

15 

6 

8 

253 

0 

4 

15 

474 

15 

7 

8 

407 

0 

4 

15 

763 

15 

8 

8 

8 

0 

4 


15 

15 

15 

9 

8 

8 

0 

4 

15 

15 

15 

10 

8 

8 

0 

4 

15 

15 

15 

11 

8 

8 

0 

4 

15 

15 

15 


14 




















15 

































































































































































n 


X2 


U2 

0 

8 

-19 

-9 

-15 

1 

424 

4109 

125 

113 

37 

168 

866 

44 

2 

141 

1296 

36 

1 

91 

657 

157 

2 

3 

11 

463 

39 

25 

15 

713 

782 

79 

4 

19 

96 

4 

23 

‘-t 

33 

131 

379 

83 

5 

45 

80 

1 

25 

83 

139 

453 

93 

6 

38 

45 

0 

4 

71 

83 

15 

7 

8 

213 

0 

4 

15 

398 

15 

8 

8 

372 

0 

4 


15 

697 

15 

9 

8 

447 

0 

4 

15 

838 

15 

10 

8 

8 

0 

4 

15 

15 

15 

11 

8 

8 

0 

4 

15 

15 

15 


8 Conclusions 

In this paper we presented a Picard iteration for solving a fixed point problem 
equivalent to a variational inequality on a cylinder. The iteration is mono- 
tonically convergent to the solution of the variational inequality with respect 
to the partial order dehned by an extended Lorentz cone. The monotone 
convergence is based on the isotonicity of the projection onto a cylinder with 
respect to the partial order dehned by the extended Lorentz cone. We plan 
to consider the following general questions in the future: 

1. Given a cone K, determine all closed convex sets C onto which the 
projection is isotone with respect to the partial order dehned by the 
cone. 

2. Given a closed convex set G, determine all cones K such that the 
projection onto C is isotone with to respect the partial order dehned 
by K. 

3. Determine the closed convex sets C for which there exists a cone K, 
such that the projection onto C is isotone with respect to the partial 
order dehned by K. 

Although the above questions are difficult to answer in general, any particular 
result about them can be important for solving complementarity problems 
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and/or variational ineqnalities by nsing a monotone convergence. Moreover, 
any snch result could be important in statistics as well, where the isotonicity 
of the projection may occur in various algorithms (see for example m- 
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